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Abstract 



We explore how nuclear modifications to the free nucleon structure func- 
tions (also known as shadowing) affect heavy quark production in collisions 
at different impact parameters. We assume that the nuclear modifications 
arise from a density dependent effect such as gluon recombination and are 
thus proportional to the local density. We calculate the dependence of charm 
and bottom quark production on impact parameter and show that density 
dependent modifications can lead to significant reductions in the heavy quark 
production cross sections in central relative to peripheral interactions. 
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Experiments have shown that the proton and neutron structure functions are modified 
by a nuclear environment. For momentum fractions x < 0.1 and 0.3 < a; < 0.7, a depletion is 
observed in the nuclear parton distributions. The low x, or shadowing, region and the larger 
X, or EMC region, is bridged by an enhancement known as antishadowing for 0.1 < a; < 0.3. 
Recently, the entire characteristic modification as a function of x has been referred to as 
shadowing. Many theoretical explanations of the nuclear effect have been proposed. In 
most, the degree of modification depends on the local nuclear density. For example, if the 
modification is due to gluon recombination , then the degree of modification should depend 
directly on the local gluon density and hence on the spatial position of the interaction within 
the nucleus. Nuclear binding and rescaling models also predict that the structure function 
depends on the local density 

Most measurements of structure functions have used lepton or neutrino beams. Typically 
these experiments are insensitive to the location of the interaction witin the nucleus so that 
the resulting structure functions are averaged over the entire nucleus. The Fermilab E745 
collaboration studied uN interactions in a bubble chamber where the dark tracks indicated 
that the interaction occurred deep within the nucleon. They found that structure function 
does vary spatially, but had little direct sensitivity to the impact parameter [Q. 

This letter discusses the effects of spatially-dependent nuclear parton distributions on 
charm and bottom production cross sections as a function of impact parameter. Assuming 
that shadowing is proportional to the local nucleon density, we study the dependence of heavy 
quark production on the nucleon structure functions and the shadowing parameterization. 

We show that the heavy quark production cross section changes significantly when this 
spatial dependence is considered. The variation with impact parameter is especially impor- 
tant because relatively peripheral collisions are often used as a baseline in searches for new 
phenomena in more central collisions 0. Large impact parameter collisions tend to probe 
the nuclear surface where shadowing is rather weak while central collisions are also sensitive 
to structure functions at the nuclear core where nuclear modifications can be large. 

The impact parameter of an event can be determined from the transverse energy pro- 
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duced in the event. The spatial dependence of shadowing on transverse energy production 
has aheady been considered 0. In Au+Au colhsions at ^/snn = 200 GeV, the center of 
mass energy per nucleon pair, a typical calorimeter can measure impact parameters as large 
as 6 = I.SRa where Ra is the nuclear radius 0. 

Heavy quark production in heavy ion collisions is dominated by gluon fusion [^^. At 
leading order (LO), charm and bottom quarks are produced in two basic processes: qq QQ 
and gg QQ with Q = c and b. The LO cross section for a nucleus A with momentum Pa 
colliding with nucleus B with momentum Pb at an impact parameter b is then 
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where Q represents the produced heavy quark. Here i and j are the interacting partons in 
the nucleus and the functions F/^ and F^^ are the number densities of gluons, light quarks 
and antiquarks in each nucleus, evaluated at momentum fraction x, momentum scale Q^, 
and location f, z in the two nuclei. The short- distance cross section, dij, is calculable as a 



perturbation series in the strong coupling constant ^^(Q^), see Ref. |T0[ . 

These leading order calculations underestimate the measured charm production cross 
sections by a constant, usually called the K factor. A similar theoretical factor, Ktt, is 
given by the ratio of the next-to-leading-order (NLO) cross section to the leading order 
result. The K factor is relatively independent of the quark pt, pair invariant mass and 
rapidity distribution, and the parton density [|lT|. However, it can vary significantly with 
quark mass and beam energy. 

Since this calculation is at leading order, we use the GRV 94 LO |jl2[ nucleon parton 



distributions, evaluated with = 1.3 GeV, = 4.75 GeV and Q = rax where 



p\ + rriQ. We also used the GRV 94 HO with the same parameters and the MRS G 



distributions M], evaluated with the same parameters for bottom and for charm, = 1.2 
GeV and Q = 2mT- The mass and scale parameters used with each set were chosen for their 
agreement with the QQ total cross section data. 



We assume that the parton densities Ff{x,Q'^,r,z) can be represented as the product 
of X and independent nuclear density distributions, position and A independent nucleon 
parton densities, and a shadowing function that contains the modification of the nuclear 
structure functions. 

Q\ r, z) = pa{s)S'{A, x, Q^ f, ^)/f (x, Q') (2) 
Ff{x, Q\h- f, z') = Pb{s')S^{B, X, Q^ 6 - f, /)/J(x, Q^) , 



where /^(x, Q^) are the nucleon parton densities, s = ^/r^~+^ and s' = y\b — r\'^ + z'"^. 
In the absence of nuclear modifications of the structure functions (no shadowing), 
S^{A, x,Q'^,f, z) = 1. The nuclear density is given by a Woods-Saxon distribution 

1 + uj{s/Ra)^ 

Pa{s) = po— ^ „ w,i (3) 

1 + exp[(s - RA)/d\ 

where electron scattering data from is used to fix the parameters Ra, d, uj and po- 
The shadowing effect is studied with two parameterizations previously used to estimate 



the effect on heavy quark production without including spatial dependence ||T5[. The first 



5*1(74, x) is based on fits to recent nuclear deep inelastic scattering data |T6|. It treats the 
quark, gluon and antiquark functions equally without evolution. The second, 5*2 (A, x, Q^) 
has separate modifications for the valence quarks, sea quarks and gluons and includes 



evolution [|T7 |. 



To include the spatial dependence of shadowing, we assume that these modifications are 
proportional to the undisturbed local nuclear rest density, Eq. (H), 

^ivs = S\A,x,Q\r,z) = 1 + iVws[5^(A,x,Q2) _ 1]^ (4) 

Po 

where A^wsis a normalization constant chosen such that (1/^4) / d^sp{s)S\i^^ = S\ At large 
radii, s ^ Ra, medium modifications weaken so that in very peripheral interactions, nucle- 
ons behave act as if they are in free space. At the center of the nucleus, the modifications 
are larger than the average value found in lepton scattering experiments. A second param- 
eterization, S*^, based on the thickness of a spherical nucleus at the collision point [0 , leads 



to a slightly larger modification in the nuclear core. Because it assumes a sharp nuclear 
surface, for s > Ra, 5*^ = 1, enhancing the shadowing at the center over that found with 
Eq. (|^). Thus for surface nucleons, Eq. predicts modifications half as strong as at the 
core while no further modifications are predicted with 5*1^. Other calculations assuming dis- 
tinct shadowing effects in different nuclear orbitals ||^ cannot be directly compared to our 
results. 

With these ingredients, Eq. (|l|) can be used to find the charm and bottom cross sections 
as a function of impact parameter. Table 1 gives the total cross sections in nucleus-nucleus 
collisions integrated over impact parameter (in units of /ib/nucleon pair) for several cases. 
As a baseline for comparison we give the results with = 1 for GRV 94 HO and MRS 
G structure functions, followed by S* = 1, 5*1, and 5*2 for GRV 94 LO. The theoretical K 
factors for each S* = 1 set are also included. We note that, with this normalization, the total 
cross section, integrated over all impact parameters, is unchanged when spatial dependence 
is included in 5*. 

Although the spatial parameterizations are normalized to reproduce the impact- 
parameter integrated cross sections, the total cross section changes with impact parameter 
in heavy ion collisions. This effect is illustrated for the GRV 94 LO distributions in Figs. 
1 and 2. To emphasize the role of the spatial dependence, the results are given relative to 
the S = 1 cross section. Figure 1 presents the charm production ratios for collisions at LHC 
(Pb+Pb at ^J]^ = 5.5 TeV), RHIC (Au+Au at ^JJ^ = 200 GeV) and the SPS (Pb+Pb 
at ^ysNN = 17.3 GeV) while Fig. 2 shows the bottom production ratios at LHC and RHIC. 

Figure 1 shows that the nuclear effect increases with the collision energy. At RHIC 
and LHC, charm production occurs predominantly at momentum fractions Xi < 0.1, in 
the shadowing region, leading to a decrease in the cross sections for both and 5*2. In 
central collisions at the LHC, the cross section decreases by a nearly 50% while at RHIC the 
decrease is about 30%. In both the impact parameter increases, somewhat stronger 

shadowing is observed in central collisions when the spatial dependence is included while for 
impact parameters larger than Ra the shadowing decreases until the ratio begins to approach 
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unity for b > 2Ra- At the much lower SPS energy, charm production typically occurs in the 
antishadowing region. Since parameterization 5*2 includes stronger gluon antishadowing than 
^i, the charm cross section is enhanced by 12% in central collisions when 5*2 is used while 
the nuclear effect with 5*1 is negligible. As the impact parameter increases, the enhancement 
decreases. Similar effects are seen with the other sets of parton distributions studied. 

We note that these cross sections are integrated over the final-state kinematics and 
that, under certain kinematic conditions, the shadowing effect can be larger. Thus, the 
rapidity and pt distributions of the produced quarks may exhibit a stronger shadowing 
spatial dependence, as shown in detail at RHIC energies in Ref. 0. E.g. when y = for 
both charm quarks, Xi = X2, providing the cleanest determination of the magnitude of the 
nuclear modification. 

Figure 2 shows the corresponding shadowing ratios for bb production at LHC and RHIC. 
The heavier b quark probes larger values of x and than the charm quark at the same 
energies and thus exhibit a different shadowing pattern. At the lower of charm produc- 
tion, the scale-dependent shadowing parameterization 5*2 produces stronger gluon shadowing 
than Si at LHC and RHIC. However, for bb production, the dependence of 5*2 reduces 
the gluon shadowing with respect to Si. Figure 1(a) shows an approximate 60% shadowing 
at 6 ~ with S2, compared with 25% for b quarks in Fig. 2(a). Meanwhile calculations 
with 5*1 show little difference between charm and bottom production. This is because at the 
typical X values at LHC energies, ~ 4.7 x 10~^ and 1.7 x 10^^ for c and b quarks respectively 



at ?/ = and pt = 0, nuclear shadowing has been seen to saturate []18| so that the larger 
X for b production has only a small effect on the Q^-independent Si parmeterization. The 
lower energy of RHIC, where x ~ 0.013 for charm and 0.048 for bottom with y = and 
Pt = 0, exhibits a larger difference between the heavy quarks, especially after evolution 
is considered. Over the kinematic range of the total cross section, bb production results in 
a 2% enhancement in central collisions due to antishadowing for S2. A 9% depletion is still 
observed for Si without evolution. Although not shown, we note that at SPS energies, b 
production is predominantly in the EMC region of x, leading again to a depletion of bottom 



production relative to that observed with no shadowing. 

Our results show that using peripheral collisions as a baseline for comparison to central 
collisions can lead to significant errors if the spatial dependence of shadowing is not taken 
into account. At RHIC and LHC, the cross section is generally higher in peripheral collisions 
than might be expected without any spatial dependence of shadowing. Interestingly this 
effect is reversed at the SPS. We note however that the spatial dependence of the shadowing 
is not evident until b > 1.2Ra, requiring studies of peripheral collisions to determine the 
strength of the effect. 

Similar impact parameter based effects should be observable in other hard processes such 
as J/ip and Drell-Yan production |T^. The Drell-Yan pairs are of special interest because 
they are produced by qq interactions at leading order, in contrast to the gluon-dominated 
heavy quark production. 

In conclusion, we have shown that introducing a very natural spatial dependence in 
nuclear shadowing changes the heavy quark production rates in heavy ion collisions, al- 
tering the expected relationship between central and peripheral collisions. This alteration 
could lead to a misinterpretation of the transverse energy dependence of certain quark-gluon 
plasma signatures due to the relationship between transverse energy and impact parameter. 
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FIG. 1. Charm quark production relative to production without shadowing, S = 1, as a function 
of impact parameter. The dashed and dotted hues show the effect with shadowing but without 
spatial dependence for Si and ^2 respectively. The spatial dependence is illustrated for Siws (open 
circles) and 52 ws (filled circles). The results are shown for (a) Pb+Pb collisions at the LHC with 
■,/sNN = 5.5 TeV, (b) Au+Au colhsions at RHIC with y^sjvjv = 200 GeV and (c) Pb+Pb colhsions 
at the CERN SPS with ^sKw = 17.3 GeV. 
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FIG. 2. Bottom quark production relative to production without shadowing, 5 = 1, as a 
function of impact parameter. The dashed and dotted hnes show the effect with shadowing but 
without spatial dependence for Si and ^2 respectively. The spatial dependence is illustrated for 
Siws (open circles) and 52ws (filled circles). The results are shown for (a) Pb+Pb collisions at 
the LHC with = 5.5 TeV and (b) Au+Au collisions at RHIC with ^sKw = 200 GeV. 
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TABLE I. Leading order cc and 66 total cross sections, in units of /xb per nucleon pair, inte- 
grated over all impact parameters, for the MRS G, GRV 94 HO and GRV 94 LO parton distribution 
as well as the theoretical K factors (Kth = ctnlo/o'lo) for S = 1. The cross sections including shad- 
owing are also given for the GRV 94 LO distributions. The quark masses and the scale parameters 
are described in the text. 
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